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Abstract 

This paper presents a complete algebraic proof of the renormalizability of the gauge invariant d = 4 operator F^^ (x) to all 
orders of perturbation theory in pure Yang-Mills gauge theory, whereby working in the Landau gauge. This renormaliza- 
tion is far from being trivial as mixing occurs with other = 4 gauge variant operators, which we identify explicitly. We 
determine the mixing matrix Z to all orders in perturbation theory by using only algebraic arguments and consequently 
we can uncover a renormalization group invariant by using the anomalous dimension matrix F derived from Z. We also 
present a future plan for calculating the mass of the lightest scalar glueball with the help of the framework we have set 
up. 



1 Introduction 

Quantum Chromodynamics (QCD) is the theory of the strong interactions which describes the force between 
quarks and gluons. At high energies, this theory is asymptotically free, while at low energies, only colorless bound 
states appear in nature due to confinement. So far, the mechanism of confinement is still poorly understood, since 
at low energies non-perturbative aspects play an important role. It is therefore of paramount importance to study 
objects which can bring us closer to the understanding of confinement. Glueballs are highly interesting candidates, 
as they are entirely composed of gluons, and therefore the gauge field itself is a crucial ingredient. For standard 
hadronic particles on the other hand, also matter fields are indispensable. Hence, glueballs have been widely inves- 
tigated, experimentally, on the lattice and in various theoretical models HI. 

So far, there is no clear experimental evidence for the existence of glueballs. If glueballs are observable particles, 
they would strongly mix with other states containing quarks. Due to this feature, a clear observation of a glueball 
state turns out to be rather difficult. However, there are already many indications for the existence of glueballs, 
and the debate is currently ongoing. It is worth mentioning here that several experiments are actually running and 
other ones are planned to start in the next future: PANDA [Tl , BES III fS] and GlueX [Tl to name only a few. 
Glueballs might also play an important role in the quark gluon plasma, a case that will be studied at e.g. the heavy 
ion collision experiment ALICE at CERN ||5l. 
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As no clear experimental data is yet available, the output of theoretical models ought to be compared with lattice 
data. In lattice gauge theories, there is no doubt about the existence of glueballs, although lattice calculations 
are still limited as they cannot determine the decay channels of glueballs. In contrast with possible experimental 
data, lattice calculations can however also consider pure gauge theory. A consensus on the lowest lying scalar 
glueball mass in the pure gauge gauge theory has already been reached : Mq++ ^1.6 GeV for SU(3) ll6l[71[8ll9l[T0l 

m. 

Many theoretical models have been investigated and compared with the lattice data. An extensive recent overview 
is given in [Ij. Historically, the first model to describe glueballs is called the MIT bag model y.2|. In this model, 
gluons are placed in a bag and confined by a boundary condition and a constant energy density B. This model, 
however, is rather phenomenological in nature. Other phenomenological models assume the gluons to have an 
effective mass U3. ,14] , which can be used to compose effective (potential) theories in which the masses of the 
different glueballs are calculated ifBlfTelflTllTSl. 

A more direct way to deal with glueballs is by identifying suitable gauge invariant operators, which carry the correct 
quantum numbers to create/annihilate particular glueball states, and then calculating the corresponding correlators 
to get information on the mass. In particular, this route is followed in the widely used QCD sumrule approach 
|fT9ll20l . For example, the operator relevant for the lightest scalar glueball is F^{x) = F^^{x), hence the study of the 
correlator (F^ {x)F'^ {y)Y One takes into account perturbative as well as non-perturbative contributions, which are 
associated with condensates and instantons 1.19. ,21j . Also in the AdS/QCD approach, glueball (correlators) have 
been investigated based on the assumption that there is an approximate dual gravity description Il22ll23l . 

In the light of such correlator studies, it would be interesting to investigate the correlator (F^ {x)F'^ [y)'^ within 
the Gribov-Zwanziger framework. The Gribov-Zwanziger action 1*241 was constructed in order to analytically 
implement the restriction to the Gribov region i2, defined as the set of field configurations fulfilling the Landau 
gauge condition and for which the Faddeev-Popov operator. 



is strictly positive, namely 

a = {a;;,8^a;; = o, 5v/''*>o}. (1.2) 

The boundary, of the region fi is called the (first) Gribov horizon. This restriction is necessary to avoid the ap- 
pearance of Gribov gauge copies in the Landau gauge f25l. Unfortunately, there are still a number of Gribov copies 
remaining, but the Gribov-Zwanziger action is so far the best approximation available. The Gribov-Zwanziger ac- 
tion is originally constructed as a non-local action. However, with the introduction of new fields, one can localize 
this action into the following form, 

-ig I d^x (^r'^Ay; + r"^A;^''; + 1 (a^2 _ ^ ^ ^ ^13) 

whereby Sym denotes the Yang-Mills action, the fields {^"i^ are a pair of complex conjugate bosonic fields 

and (ro^'^j 0)"') are anticommuting fields needed to localize the original non-local action. The parameter y is fixed 
by a gap equation and implements the restriction to the Gribov region. This y introduces a mass scale into the 
theory which can consequently give rise to a nonvanishing pole in the glueball correlator. By taking into account 
the dynamics of the fields ((p,(p,co,co), it became clear that a second mass scale emerges quite naturally, and 
this can also enter the glueball correlator expression. For more details on the Gribov-Zwanziger action, its 
renormalization and the dynamics of its constituent fields, we refer to | 



Once the operator F^y is introduced, the issues of renormalization and mixing complicate matters at the quan- 
tum level. If one wants to investigate F^^ in a renormaUzable setting, one should introduce this operator into the 
action by coupling it to a source q, and then renormalize that action. Therefore, before scrutinizing the more com- 
plicated Gribov-Zwanziger case, it is instructive to first completely investigate F^y within the usual Yang-Mills 
theory, quantized in the Landau gauge. Moreover, as F^y is not a renormalization group invariant, one could look 
for a renormahzation group invariant operator, containing F^^, since renormalization group invariance is beneficial 
when looking at physical quantities, in casu the glueball mass. This renormalization group invariant will turn out 
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to coincide with the trace anomaly. However, we would like to avoid a direct use of the trace anomaly, as the 
renormalization of the trace anomaly itself, through that of the energy momentum tensor, is rather difficult and 
sometimes tricky f28l. Therefore, we shall focus on the direct renormalization of F^^ in the Landau gauge, which 
also turns out to be far from trivial as a mixing with other (non gauge invariant) operators occurs. In 1974, [291 
described the first attempt towards the renormalization of gauge invariant operators. In this paper, the renormal- 
ization of f^y at zero momentum was investigated to the first loop order, and the renormalization group invariant 
containing F^^ determined. However, the paper |29 1 focused only on the integrated operatoiQ / A'^xF^^ . In addition, 
a generalization to the more complicated Gribov-Zwanziger case does not seem straightforward to implement in 
the language of |29|. Also, no clear proof of the higher order renormalization of the nonintegrated operator Fp^{x) 
can be found. Notice that the passing from the integrated to the nonintegrated operator is not trivial, see |30|, §12.6, 
and references in ||3T1 . 

In ll32l[33l[34l . one has elaborated on the structure of the mixing matrix for the more general case of non-integrated 
gauge invariant operators, while in fSFl a simplified proof of the renormalization of gauge invariant operators 
has been given from the perspective of the BRST cohomology of Yang-Mills gauge theories. In the light of this 
reference ll35l . the last paper published on these issues focusses in particular on some infrared subtleties 1361 . One 
can appreciate the intrinsic difficulties arising when studying gauge invariant operators at the quantum level by 
noticing that in [i37i . results of 1,32] were used, while one of the authors of [371 quotes the same paper |32 | again 
in later years in ll36l . mentioning that he finds the proof of ll32l "very hard to understand". 

Although the cohomological proof of ||35l is of full generality, the results are of an abstract nature. In an oversim- 
plifying nutshell, it was shown that each BRST invariant operator at the quantum level with ghost number zero, 
can be written as a strict gauge invariant operato:[l plus BRST exact piece, modulo terms that vanish upon using 
the equation of motions. For practical computations, it is however not sufficient to know which type of operators 
occur in the renormalization process, but also the explicit knowledge of all these operators is necessary. 

For completeness, we also mention (3W\, concerning the renormalization of F^^ in the background gauge formaUsm, 
which would however be of little use when looking at the Gribov-Zwanziger generalization. 

Based on all the foregoing arguments, we have found it instructive to present in this paper a clean analysis of the 
renormalization of the operator F^^ to all orders of perturbation theory. The proof shall be given in the framework 
of algebraic renormalization [39], and we shall retrieve a renormalization matrix, restricted by various Ward iden- 
tities. Next to suitable adaptations of the usual Landau gauge Ward identities, we also identify a new powerful 
identity, relevant in the discussion of the renormalization matrix, which form shall be in perfect agreement with the 
argumentation given in (30 36]. We stress that our analysis is purely algebraic, and in this sense differs from the 
argumentation given in [30, 36 1. Moreover, we shall also be able to completely fix the mixing matrix to all orders 
and this without calculating any loop diagram. Multiple checks will be presented, which will confirm the results. 
For example, we shall recover in an independent fashion well-known nonrenormalization relations in the Landau 
gauge [39], here stemming from the renormalization analysis of F^^. 

In summary, in section II we shall give an overview of the 3 different classes of operators which can mix with F^y 
before going into the detailed algebraic renormalization of F^^. In section III, the mixing matrix will be determined 
to all orders in perturbation theory and, armed with this result, we shall be able to construct a renormaUzation 
group invariant in section IV. We end this paper with a discussion in section V. 



2 Renormalization of the Yang-Mills action with inclusion of the operator F^^ 
2.1 Introduction 

The most natural way to study the lightest scalar glueball is by determining the correlator (^F^ {x)F^ {y)y This 
correlator can be obtained by adding the operator F^^ to the ordinary Yang-Mills action by coupling it to a source 
q{x). 

Indeed, the action we start from reads, 

En.. = jd'x\F^, + jA'x{b"d,Al + c^%Dfc')+jd\'iF^,, (2.1) 

" V ' 

5gf 




1. Or the operator at zero momentum. 

2. This means containing only the field strength and covariant derivatives. 
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whereby Sgf is the Landau gauge fixing part. For the benefit of the reader, let us already mention that this action is 
BRST invariant, 

*2n.r. = 0, (2.2) 

with all the BRST transformations of the fields and the source given by 

sc" = b", sb" = 0, sq = 0, (2.3) 

and s nilpotent. 



= 0. (2.4) 



In this fashion, the correlator is given by 

5 5 



1=0 



{FHx)FHy)) , (2.5) 



with Z*^ the generator of connected Green functions. However, it will turn out that the action (12.11 ) is not renormal- 
izable. Indeed, as the operator F^^ has mass dimension 4, it could mix with other operators of the same dimension. 
The question arises which kind of extra operators we need to consider. 

2.2 3 classes of operators 

In general, we can distinguish between 3 different classes of dimension 4 operators. Firstly, the class Ci contains 
all the truly gauge invariant operators. These are the BRST closed but not exact operators like F^^. These are 
constructed from the field strength F^^ and the covariant derivative D^''. Secondly, the class C2 consists of BRST 
exact operators, e.g. s{c"dfjA'p. The third class C3 contains operators which will vanish upon using the equations 

of motion, e.g. A^-gp, with S = Sym + Sg[. 



Now, one can intuitively easily understand that these 3 different classes will mix in a certain way H30II36I . Firstly, 
bare operators from the class C2 cannot receive contributions from gauge invariant operators (Ci). Indeed, taking 
the matrix element of a bare BRST exact operator from C2 between physical states will give a vanishing result, if 
there would be a renormalized gauge invariant contribution from Ci in its expansion, there would be a nonvanishing 
contribution, clearly a contradiction. Secondly, as a C3 operator will vanish upon using the equations of motion, 
while a Ci- and a C2 operator in general do not, a C3 operator cannot receive corrections from the Ci and/or C2 
class. 

Thus, the mixing matrix will have an upper triangular form. 





Z^^ Z^^ \ \ L \ , (2.6) 
whereby f , L, J{ are operators belonging, respectively, to the C\, C2 and C3 class. 

We shall however not use these observations, and only rely on a formal algebraic analysis ll39l . All constraints on 
e.g. the mixing matrix should be encoded in the Ward identities. 

For further use, let us elaborate a bit more on the equation of motion like terms, using a scalar field for nota- 
tional simplicity. A term ^ |^ shall give rise to contact terms when taking expectation values. Using partial path 
integration, one finds 

(p(xi)(p(.X2)...(p(x„+i)^3;373 ) = / [d(j)](p(xi)(p(x2)...(p(x„+i)^3;373e-^ 

e-^ 



- j [d*](p(-*:i)(p(-*:2)...9(-««+i)g^ 
f 6 

1 ['**]5;^['P(-^i)^(-^2)---^(-^«+i)]^ 



89(3/) 

n+l 

^ h{xk ~ y) ((p(xi )(p(x2) . . .v?{xk-i)v?{x,,+ i) . . . (p(x„+i )) . (2.7) 



k=\ 



4 



We used the symbolic notation /[d(|)] for the integration over all the present fields. Introducing the Z-factors for 
the fields cp , one also learns that (f{y) does not need any renormalization factor, and thus that it is finite when 

introduced into correlator^ Moreover, if xi^ =/= y, k = the l.h.s. of ( 12. 7t will vanish as the r.h.s. does. On 

the other hand, it is easily recognized from ( 12.7b that the integrated operator / d'*y(p(y) is nothing more than a 
counting operator when inserted into a correlator, i.e. 

Atp(3')g^y = «((p(xi)(p(x2)...(p(jc„)) . (2.8) 
2.3 The starting action 

We can now propose a more complete starting action than (12.1b . Besides the gauge invariant operator F^^ belonging 
to the first class Ci, we also introduce the BRST closed operator s{dcA) = i(3^c"A^), coupled to a new dimen- 
sionless source rj. As we want this new source to only enter the cohomological trivial part of the action, we shall 
introduce a BRST doublet (A,,r|), 

sTi = A., (2.9) 

and add the following term to the action ( 12.1b . 

j A^sivi^d^l,) = |d4jc(Hc«A;;+ri(a^/7«A;;+a^c"z)fc*)). (2.10) 

The BRST doublet-structure is highly useful in order construct the most general invariant counterterm lf39ll . 
Hence, the classical starting action is given by 

5ci = 5ym + 1 A^x (b"d,Al + rd.Dfc'') + j d^xq ^F^^ + 1 d^xXd.c^'A; + J d^x^ (d^fA"^ + d,c"Dfc') . 

(2.11) 

Later in this paper, we shall also introduce the equation of motion terms from class C3. Notice that in principle, also 
i(c''3^Ap is another independent d — 'X BRST exact operator which could play a role. It shall however turn out that 
the renormaUzation analysis closes without this operator, therefore we decided to immediately discard it. 



We can now proceed with the study of this action, using the formalism of algebraic renormalization 11391 . A first 
step is to introduce a term 5ext. 



I d'x {-K; (p,cr + \gL"f^'c'A , (2.12) 



'^ext — 

needed to define the nonlinear BRST transformations of the gauge and ghost fields. K'^ and L" are two new sources, 
invariant under the BRST symmetry s. Therefore, the enlarged action is given by 

E = Sym + I d^x {b"d,A; + T^d.ofc'') + J d^x (^-K"^ {D,c)" + ^-gL''f'"'c"c'^ 

+ j d^xq^-F^^ + j d4xHc"A; + jd^'xn (d^fA^^ + d^T^'Dfc') , (2.13) 

and it will reduce itself to equation (12. lib , once the sources K^^ and L" are set to zero at the end. Likewise, also X 
can be set to zero at that point. 

A second step in the process of algebraic renormalization is to determine all the Ward identities obeyed by the 
action (12. 13b . which we have summarized here: 

• The Slavnov-Tayloridenitity: 



3. The implied limit x,i^[ y might seem problematic due to the appearance of a S(0) in the last term of the r.h.s. of 12. 7t . However, 5(0) = 
in dimensional regularization. 
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The Landau gauge condition: 



^ = 9,A«-a,(TiA«). (2.15) 



The modified antighost equation: 



g+a^-O.*^ = 0. (2.,6, 



• The ghost Ward identity: 

/d^-(5^+^r'^(^5^))^ = ^/dV^'«^;,-^V). (2.17) 

The term A"j, being linear in the quantum fields A^, c", is a classical breaking. 

• The extra integrated Ward identity: 




(2.18) 



Apart from some small adaptations, the first 5 symmetries are similar to the ones in the ordinary Yang-Mills action. 
Moreover, we also find an extra Ward identity w.r.t. the new doublet (A,,r|). This last identity will enable us to take 
into account in a purely algebraic way the effects related to the composite operators coupled to the sources (X,ri). 
We underline here that this is the power of the algebraic formalism: by a well chosen set of sources to introduce the 
relevant operators, one can hope to find additional Ward identities which, in turn, will constrain the theory at the 
quantum level, including the characterization of the most general counterterm. As such, a good choice of sources 
can considerably simpUfy the renormalization analysis. 



2.4 The counterterm 

When we turn to the quantum level, we can use these symmetries to characterize the most general allowed invari- 
ant counterterm I,'^. Following the algebraic renormalization procedure, is an integrated local polynomial in the 
fields and sources with dimension bounded by four, and with vanishing ghost number. The previous, nonanoma- 
lous. Ward identities imply the following constraints on 1/^: 

• The linearized Slavnov-Taylor identity: 

!BeZ'' = 0, !b| = 0, (2.19) 



J y5/:«5A« 



, ^4 , 5E 8 5E 8 8E 8 8E 8 8 , 8 \ 



The Landau gauge condition: 



The modified antighost equation: 



5E^ 



0. (2.21) 



S+S.IK-O.n)^ = 0. (2.22) 



The ghost Ward identity: 



8 +gfabc(-h^ ) jj^c ^ o_ (2.23) 
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A" 


c" 


c" 


b" 


dimension 


1 





2 


2 


ghost number 





1 


-1 






Table 1: Quantum numbers of the fields. 





k; 


L" 


q 


11 


X 


dimension 


3 


4 











ghost number 


-1 


-2 








1 



Table 2: Quantum numbers of the sources. 



The extra integrated Ward identity: 



To construct the most general counterterm. Tables [T] and |2] listing the dimension and ghost number of the various 
fields and sources, are useful. There is however one subtlety concerning counterterms quadratic (or higher) in 
the sources. Only looking at the dimensionality, the ghost number and the constraints on the counterterm, it is a 
priori not forbidden to consider terms of the the form {q^ . . .), (r|^ . . .), (^r| . . .), {q^ . . .), etc., i.e. terms of quadratic 
and higher order in the sources. If these terms are allowed, an infinite tower of counterterms would be generated 
and it would be impossible to prove the renormalizability of the action as new divergences are being generated, 
which cannot be absorbed in terms aheady present in the starting action. However, we can give a simple argument 
why we may omit this class of terms. Assume that we would also introduce the following term of order q^ in the 
action. 



/ 



dV^- (2.25) 



Subsequently, when calculating the correlator, this term would give rise to an extra contact term contribution, 

[d(j)]e"^ 



5q{z) 8q{y) 



= {F\z)FHy)) +?>{y^z){F\y)) . (2.26) 

q=0 " V ' V ' 

term due to part in q term due to part in 



As eventually we are only interested in the coiTelator for z 7^ y, we can thus neglect the term ( 12.25b . In fact, when 
looking at the case z — y,we should also couple a source to the novel composite operator F'^, which is not our 
current interest. We can repeat this kind of argument for all other terms of higher order in the sources. 

There is one exception to the previous remark: we cannot neglect higher order terms of the type {Kq. . .) and 
{Kr[ ...) due to the modified antighost equation. 

The second term of this equation differentiates the counterterm w.rt. the source K'^, while the first term w.r.t. the 
field c". Therefore, for the construction of the counterterm fulfilling all the constraints, we still need to include 
terms of order Kq and Kr[, as when deriving w.r.t K"i, these terms will become of first order in the sources, just 

as the term 0= However, at the end, after having completely characterized the counterterm, we can ignore this 
class of terms again. 

We are now ready to construct the counterterm. Firstly, making use of general results on the cohomology of gauge 
theoiies [ 39 1 , the most general integrated polynomial of dimension 4 in the fields and sources, with vanishing ghost 
number and which takes into account the previous remarks on the terms quadratic in the sources, can be written 
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+ «r I d^x f^biq{K; + d,r )Al + b2qc^d,A"^ + b^qb^c^' + fe4?g/"''c"c^c' j 

+«i / d^^Jc j di A.c"c" 1 . (2.28) 



Secondly, we can simplify this counterterm by imposing all the constraints ( |2.19t - (l2.24l i. After a certain amount 
of algebra, we eventually obtain 

r- = a,jd\\F^, + b,fA'xlF^,+a,fA'xi^l^-^+A'^^ 

+bi I d'^xqi A^4^ + K';d,c" + d,rd,c" I , (2.29) 



with 



Sym = ^J d'^xqF^^. (2.30) 

Now that we have constructed the most general counterterm obeying all the Ward identities, we can neglect, as 
previously described, the term in Kq. Therefore, the final counterterm becomes. 



ao I d'x^-F^,+bo I d'x^F^^+a, j d^x ( +a;;^ +/r;3,c" +a,c«a,c« -iiS.^'V" 



+bx / d'xq I Al^+d,c^V . (2.31) 



2.5 Introducing the equations of motion 

We still have to introduce the equations of motion as described in Section IZ21 as these can enter the operator f . 
So far, we have found an action E with corresponding counterterm E' . Let us perform the linear shift on the gluon 
field A", 

A;^A; + 7A;, (2.32) 

with J{x) a novel local source. This way of introducing the relevant gluon equation of motion operator shall 
turn out to be very efficient, as it allows us to uncover the finiteness of this kind of operator Indeed, this shift 
basically corresponds to a redefinition of the gluon field, and has to be consistently done in the starting action and 
counterterm. Performing the shift in the action gives rise to the following shifted action E', 

= 5ym + 1 d^x {b"d,A; + c"d,Dfc') + 1 d'x (^~k; {D,cf + ^-gL^f'c'c'^ 

+ j d\q^F^^ + j d\lJd,rAl + j d^xil {d,b"A"^ + d,c^Dfc'') 

+ j dVA«^+ j A''xj[-d,b"Al + gUtA'lc'^d,:c"] , (2.33) 

3{ 
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where we see the relevant gluon equation of motion term, 9i , emerging. Again, we have neglected higher order 
terms in the sources, as the argument (|2.26t still holds. Analogously, we find a shifted counterterm, 

+a 1 1 A^xJ (2A^8^8vA^ - IA'^/A^ + 9g/„,,A;:A^a^A^ + 4g V«/-c/.rf.A;;A^A^A^) , (2.34) 
where one can neglect again the higher order terms in the sources. 

One could also introduce the other similar equation of motion terms, by introducing linear shifts for the h" , c", c" 
fields. However, the corresponding equation of motion operators will not mix with and are therefore unneces- 
sary to establish the renormalizability of the action (12.331 1. 



2.6 Stability and tlie renormalization (mixing) matrix 

Finally, it remains to discuss the stability of the classical action, i.e. to check whether E" can be reabsorbed in the 
classical action E' by means of a multiplicative renormalization of the coupling constant^, the fields {(|) = A,c,c,i>} 
and the sources {<!> = L,K,q,T\,X,J}, namely 

l![gA^)+hT^' = ngoM,^o) + 0{h^) , (2.35) 
with h the infinitesimal perturbation parameter. The bare fields, sources and parameters are defined as 

_ ■y^l'^ a 

(2.36) 

We also propose the following mixing matrix, 

/ '^qq ^qT\ ^qJ \ / ^ \ 

(2.37) 



which will represent the mixing of the operators [F , L and y{ . If we try to absorb the counterterm into the original 
action, we ultimately find, 

Zg = 1 — H — , 







ja 


= ZiL\ 




= 'Zqq-, 


110 




^0 


= Zy^, 








2 



Zf = l+h 



(y+fli), (2.38) 



and 



'A ' 



Zk — z}J^ , 

Zl = Z'f, (2.39) 
results which are known from the renormalization of the original Yang-Mills action in the Landau gauge |[39ll . 
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In addition, we also find the following mixing matrix 



T^q 



'--qK\ 



^qj 



Zy^ Zjy 

and for completeness, the Z-factor of X reads. 



hbi 
hbi 




as the counterterm does not contain the source A,. 



7-1/2-1/2 



(2.40) 



(2.41) 



Once having this mixing matrix at our disposal, we can of course pass to the corresponding bare operators. For 
this, we shall need the inverse of the mixing matrix ( |2.40l l. 



(2.42) 



Now we can determine the corresponding mixing matrix for the operators, since insertions of an operator corre- 
spond to derivatives w.r.t. to the appropriate source of the generating functional Z'^{q,r\,J). In particular. 




To = 
Jo = 



5Z'(^,r|,y) 5^ 5Z^(^,ri,7) 5ri 5Z'(^,ri,y) 5/ 5Z'(^,r|,7) 



5^0 



5^0 5^ 



5^0 5-/ 



7 7 7 



(2.43) 



and similarly for and 9iQ. In summary, we find 



To \ 




/ '^qq 


^Jq^qq 


^Jq^qq 






I 


1 





^0 / 









1 




(2.44) 



From this matrix we can make several interesting observations. Firstly, we see that the operator \F^y {— f) indeed 
required the presence of the BRST exact operator L and of the gluon equation of motion operator y{ as these 
operators are "hidden" in the bare operator f . Secondly, we do retrieve an upper triangular matrix, in agreement 
with the earlier description in ( I2.6l l. Moreover, we also find that the BRST exact operator L does not mix with ^ , 
a mixing which is in principle allowed, but has a Z-factor equal to 1 . This can be nicely understood: the integrated 
BRST exact operator is in fact proportional to a sum of two (integrated) equations of motion terms. 



/ 



- / d\x(b"duA' 



-rd,Dfc') = - 



' 8b" 5c" 



(2.45) 



and therefore it does not mix with other operators, just like y{ . 



3 The mixing matrix to all orders 

In this section, we shall demonstrate that we can determine the mixing matrix (12.441 1 exactly, i.e. to all orders of 
perturbation theory. For this purpose, we shall follow the lines of ||40]| . suitably adapted to the gauge theory under 
study. We start with the following most general (n + 2m + r) -point function defined as, 

g"^^'"'^''{xu. . . ,Xn,yi, . . . ,yn,yu - ■ ■ ,yn,zu - ■ ■ ,Zn) = {A{xi) . . .A{x„)c{yi) . . .c{y,„)c{yi)c{y,n)b{zi) . . .b{zr)) 

= J [d(if]A{xi) . . A{x„)c{yi) . . .c{y„)^yi)c{y,n)bizi) . . .bizr)e-^ , (3.1) 

with the action S given by 

S = ^YM+'Sgf (3.2) 
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We have immediately assumed that there is an equal amount of ghost and antighost fields present as in any other 
case, the Green function ( 13.11 ) would be zero, due to ghost number symmetry. Subsequently, from the definition 
( 13. Il l, we can immediately write down the connection between the renormalized Green function and the bare Green 
function, 

From the previous equation, we shall be able to fix all the matrix elements of expression ( I2.44l i. based on the 
knowledge that 

A^n+2m+r 

^3.4) 

dg2 

is finite. 

We start by applying the chain rule when deriving the right hand side of equation (13.3b w.r.t. g^. We find, 

- ^-^^A A- ^b +—^^c H ^0 +Za -^^b ^0 A- "a^yo , 

(3.5) 

where we have replaced (n + 2m + r) with n' as a shorthand. Next, we have to calculate all the derivatives 
w.r.t. g^. 

• Calculation of ^ 

ag- 
in dimensional regularization, with c/ = 4 — e, one can write down 

gl - A'%V- (3.6) 



Hence, if we derive this equation w.r.t. g^, 

dgl ,dZl . ^ , , / ainZ2 1 , 



From the previous equation, we still have to determine which can be extracted from equation (I3.6l l. 

Deriving this equation w.r.t. /j gives, 

M = e^^zy+/^^^-§^g'+^^z'M=0, (3.8) 

were we have applied the chain rule again. We can rewrite this equation making use of the definition of the 
P-funtion 

= -e^' + P(/), (3.9) 
where we have immediately extracted the part in e, and we obtain, 

'^""^ - ' l] -\(—^)- (3.10) 



(3.11) 



If we insert this result into expression ( 13.7b . we ultimately find 

M _ -egl 



dg^ -^g2 + ^{g 



2\ ■ 



-n/2 



Calculation of 

dg- 

The next derivative w.rt. g^ can be calculated in a similar way. We start by applying the chain rule, 

_ _„V^<! _ _,^^_„/23lnZ'/^ 

^A 



dg2 - "^1/2 3^2 - " 3^2 • (3.12) 
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Next, we derive from the definition of the gluon anomalous dimension, 



Ya = A* 



3lnZy^ _ dg^ dXnZ^f 
d/j ^ d/j dg^ 



dlnZ 



1/2 



dg 



2 



From expression ( 13.121 ) and ( 13.131 ), it now follows 



(3.13) 



^^A ^ _ ^-"/2 Ya 

dg' -^8' + m 



(3.14) 



^7— /J! 

Calculation of 

Completely analogously, we find with the help of the anomalous dimension of the ghost field, 

dz'J' 1/2 



(3.15) 



(3.16) 



Calculation of — 



Finally, from 



dz, 



1/2 



„l/2 



(3.17) 



we deduce 



dz, 



-r/2 



dg' 



rZ 



-n/2 



Ya 



(3.18) 



Taking all the previous results into account, we can rewrite expression ( 13.51 ) 
dg"' 



Ag' 



_8_ 



e^o T-2 - "Ya - 2mY,. - ry/, 



(3.19) 



The right hand side of ( |3.19l l still contains bare quantities, which we have to rewrite in terms of renormalized 
quantities. Notice also that we would like to get rid of the field numbers n, m and r as the mixing matrix will 
evidently will be independent of these numbers as they are arbitrary. 

We shall now alter the right hand side of equation ( |3.19t by calculating and by removing the fields numbers. 

After a little bit of algebra, we obtain. 



8(e 



dgl 



2gl 



Ao{y) 



5S 



5AoW 



2gl 



My)dAo{y)] e 



(3.20) 



Consequently, deriving the n point Green function g'^ w.rt. gQ will result in several insertions in this Green 
function. 



dgo 



1 



4 r2^«i^°(^)5AoW 



55 



(3.21) 



where we have introduced a shorthand notation, e.g. 



^^0 



n'lF^iy) 



FoHy) 



A{xi).. .A{x„)c{yi) . . . c(y„)c(yi )c(y„)fe(zi ) . . .b{zr) 



(3.22) 
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The field numbers can be rewritten by inserting the corresponding counting operator. If we start by inserting the 
counting operator for the gluon fields n, we find 



5S 



5AoW 



(3.23) 



as derived in equation (12. St . We can derive analogous relations for the other counting operators. 



55 1 




§co(y) J 


■ = m^o 






5^0 (y) J 





(3.24) 



Taking all these results together, expression ( 13.191 1 now becomes. 



•7 '^I'^-y—m-y '"/^ 



(3.25) 



(Jo 

We have again introduced a notational shorthand for the equation of motion operators, with k = bo-§r- and 



Kq = cq-^ and with f and already defined before. 

In the last part of the manipulation of the n'-point Green function we reexpress all the operators again in terms of 
their renormalized counterparts, thereby writing all the divergences explicitly in terms of 8. Firstly, we can reabsorb 
the Z-factors into (^q to find. 



^—^ I -e (g"' {My)} - {^oW} + \g"^ {loiy)} 



-lA{g^)g"\^,{y)]-2lc{g')g"\%i^{y)}-it{g')g"\io{y)} 

Secondly, we parametrize the mixing matrix ( I2.44l i. 



(3.26) 



^0 









1 



L 
9{ 



(3.27) 



Here we have displayed the fact that the entries associated with a{g^,t) and b{g'^,t), which represent a formal 
power series in g^, must at least have a simple pole in e. We recall that the integrated operator Lq is proportional to 
the sum of the two counting operators / /o and / see expression ( 12.451 ). Therefore J lo~ J l and j "Ka — J Ki- 
Inserting all this information into expression ( 13.261 1 yields. 



dg" 



1 



j d^ g-'^^{-^-a)J{y)-bI{y)-b^{y)+bM{y) + ^-M{y)-^IiJ) 



{g')9{{y)-2Ug')x{y)-yt{g'')i{y)} 



(3.28) 



giving us the final result from which we shall be able to fix the matrix elements of expression ( I2.44l i. 

As the left hand side of our final expression ( 13.28b is finite, the right hand side is finite too. Therefore, the following 
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coefficients must be finite, 

7 
I 



-Z-a 1 (1+a/e) 



-e/2-/7-Yi(^2) 1 l+2(/7 + y,,(^2))/g 



£/2 + fc-YA(g') ^ 1 l+2(/7-YA(g'))/e 

-eg2 + p(^2) 2g2 l-P(g2)/(eg2) ■ 



K : -'-'Jfl (3.29) 
seen as a power series in g^. This can only be true if 

%',e) = YA(g')-^^--Y/.(g')-^^--2Yc(g')- (3.30) 

2 §2 2 g2 

In fact, this last equation reveals a connection between the anomalous dimension of A and b, and between the 
anomalous dimension of A, g and c, namely 



lA+lb = 0, 



JA +2yc = (3.31) 



1/2 

These relations are well-known to hold in the Landau gauge, since Z^Zf, = 1 and Z^Z^ '^g = 1' derived from 
the algebraic renormalization analysis, which leads to equations (|2.38l l and (|2.39t . 

In summary, we have completely fixed the mixing matrix in term of the elementary renormalization group func- 
tions, and this to all orders of perturbation theory, 

' 1 I I i I . (3.32) 





^0 z 

In addition, as a check of this result, we have also uncovered two relations, ( 13.31b . between anomalous dimensions 
which must hold for consistency. These correspond to ( 12.381 1 and ( |2.39l l. which are well-known nonrenormalization 
theorems in the Landau gauge. 



4 Constructing a renormalization group invariant 



As a last step, we can now look for a renormalization group invariant operator by determining the anomalous 
dimension Y coming from the mixing matrix Z. We define the anomalous dimension matrix Y as 



For the calculation of F, we require 



A* 



so we obtain. 



1 - 



.2 9(P/r) 






= ZF. 



-i-^g' 



2 dye 






(4.1) 



= --(-e^2^p(^2)) 



mg^ 
dg^ 




(4.2) 



(4.3) 
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which is indeed finite, a nice consistency check. This matrix is then related to the anomalous dimension of the 
operators: 

dz dx 

Xo ^ ZX =^ = fj—X+Z/j— 
a/j a/j 

dX 

^ ,j— = -FX. (4.4) 

We now have all the ingredients at our disposal to determine a renormalization group invariant operator. We are 
looking for a linear combination of f , L and y{ which does not run, 

Ij^lkf +£g +m:N:] = 0, (4.5) 

whereby k, i and m are to be understood as functions of g^. Invoking the chain rule gives 

dk 7d(Q/g^) .djc ydjc de dm 

^—f-kg^^^^T+kg^^L+kg^^9{+p—L+iJ—^ = 0. (4.6) 
OjU ag^ ag^ ag^ d/j ofi 

This previous equation results in two differential equations, 

^£ = m, 
which can be solved by, 

\iig^)=mig^)^-y,ig'). 

In summary, we have determined a renormalization group invariant scalar operator ^ containing 5" . Explic- 
itly, 



4 g2 



-F^v~yc{g') {A;d,b" + 3,0^0^0") -y,.{g^)A;—, (4.7) 



without having calculated any loop diagram. Moreover, this invariant is equal to the trace anomaly which is 
expected as is also a d = 4 renormalization group invariant 1^37 j . 



5 Conclusion 

In this paper, we have provided a detailed analysis of the renormalization of the non-integrated operator F^y in 
Yang-Mills gauge theory in the Landau gauge, and this to all orders. We have shown that this operator mixes 
with two other operators: the BRST exact operator s{c"dfjA'p and a gluon equation of motion operator, A'^^- We 
have composed the corresponding renormalization matrix Z and we have been able to determine this matrix to 
all orders in perturbation theory. Several checks have confirmed the results. We have recovered the well known 
non-renormalization theorems in the Landau gauge, relating the gluon, ghost, auxiliary field and coupling constant 
anomalous dimension. In addition, we have calculated the anomalous dimension matrix using Z, which was nicely 
finite, and we have been able to construct a renormalization group invariant containing F^y. 

When turning to physical states, the BRST exact term and the equation of motion term will drop and be no longer 
relevant. However, we have paved the way for more complicated actions as this framework is very solid. For 
example, in the Gribov-Zwanziger action, when investigating the operator F^y, we also expect a mixing. Again, 
we would expect three different classes of operators, the gauge invariant, the BRST exact and the ones which 
vanish upon using the equations of motion. However, the Gribov-Zwanziger action is no longer BRST invariant 
Il27ll4n . Fortunately, the framework we have set up can be saved as one can embed the Gribov-Zwanziger action 
into a "larger" BRST invariant action which will reduce to the known Gribov-Zwanziger action in the so-called 
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physical limit ||27l . The embedding into a BRST invariant action seems to be the crucial tool to correctly identify 
the relevant, renormalizable, gauge invariant operators in the Gribov-Zwanziger case, as the physical limit of their 
more general counterparts constructed from the larger action. In this case, the BRST exact operators shall no longer 
be automatically irrelevant when taking the physical limit, as the BRST symmetry is softly broken by the restriction 
to the Gribov region. Therefore, we expect a non-trivial mixing to occur when we transfer to the physical reality of 
glueballs, i.e. when calculating the corresponding correlator Hence, it will certainly be of the utmost importance 
to take all the possible mixings into account in order to be able to construct a renormalization group invaiiant 

ma. 

We did not include fermions in our analysis. The cohomological treatment of massless fermions is not that com- 
plicated, see (39]. However, additional equation of motion terms, as well as the derivative of the singlet vector 
current can be relevant. In addition, including non-degenerate fermion masses would complicate matters, due to 
the introduction of a fermion mass mixing matrix. As we already mentioned, our interest is and will be mainly 
focused on pure gauge theories ll42l . 
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